NASA TM X-73635 


Vh *. 


* 


NASA TECHNICAL 

MEMORANDUM 


NASA TM X-73635 


r IBASA-TM- X-73635) VELOCITY AND TENPEBATUBE 
PROFILES IN NEAB-CBITICAL NITBOGEN FLOBIN3 
PAST A HOBIZCNTAL PLAT PLATE (NASA) U p HC 
A02/HF A01 CSCL 20a 


N77-30999 

Unclas 

42006 


{ 

\ 

VELOCITY AND TEMPERATURE PROFILES IN NEAR-CRITICAL 
NITROGEN FLOWING PAST A HORIZONTAL FLAT PLATE 

by Robert J. Simoneau 
Lewis Research Center 
Cleveland, Ohio 44135 

TECHNICAL PAPER to be presented at the 
Seventeenth National Heat Transfer Conference 
sponsored by the American Society of Mechanical Engineers 
Salt Lake City, Utah, August 14-17, 1977 


f. v V 


ABSTRACT 


Subscripts 


Boundary layer velocity and temperature profiles 
were measured for nitrogen near its thermodynamic 
critical point flowing past a horizontal flat plate. 
The heated surface was oriented both facing upward 
and downward. The results were compared to earlier 
work in which measurements were made for vertically 
upward flow. The boundary layer temperatures ranged 
from below to above the thermodynamic critical tem- 
perature. For wall temperatures below the thermo- 
dynamic critical temperature there was little varia- 
tion between the velocity and temperature profiles in 
the three orientations. In all three orientations 
the point of crossing into the critical temperature 
region is marked by a significant flattening of the 
velocity and temperature profiles and also a decrease 
in heat transfer coefficient. As the heat flux and, 
consequently, wall temperature are further Increased 
significant changes occur In the velocity and temper- 
ature profiles. Examination of ncar-critlcal heat 
transfer in these three flow orientations offers in- 
sights into the relative role of buoyancy forces in 
this regime. 

NOMENCLATURE 

C p specific heat capacity, J/g K 
g acceleration of gravity, cm/sec 2 

Gr x Graehof number, p^Cp^ - P w )x 3 g/n* 

h heat transfer coefficient, W/cm 2 K 
Nu x Nusselt number, hx/X 

P pressure, MPa 

Pr Prandtl number, nC p /X 

q heat flux, W/cm 2 

Re x Reynoldc number, p.u^x/n* 

St x Stanton number, h/UcP^Cp^ 

T temperature, K 

U velocity, m/sec 

x distance from leading edge of plate, cm 
y normal distance from the plate, mm 

X thermal conductivity, W/cm K 

n dynamic viscosity, g/cm-sec 

p density, g/cm 3 


c critical point conditions 

calc calculated 

exp experimental 

w wall conditions 

x conditions at axial position x along the 

plate 

® free stream conditions 

INTRODUCTION 

Because of low research activity, designers are 
reluctant to work with fluids near their thermodynamic 
critical point; however they remain of interest in de- 
sign. The space Industry stores its cryogenic liquids 
above their critical pressure to avoid having -apor 
present In a low gravity field. Supercritical power 
cycles have some attractive features and the cooling 
of cryogenic equipment, such as cryomagnets, is a con- 
tinuing need. Much of the low confidence for design 
analyses stems from a lack of detailed knowledge of 
the mechanisms at work, especially in forced convec- 
tion. Surveys by Hendricks, et al. (1), Hall (2), and 
Petukhov (3) yield a fairly comprehensive picture of 
the status of near-critical heat transfer. Experi- 
ments probing the details of forced convection in near- 
critical fluids have been limited. An early paper by 
Wood and Smith (A) revealed the M-shaped velocity pro- 
file which obtained its maximum between the wall and 
tube centerline. Sabersky and Hauptmann (5) added 
some detail on flow over a flat plate, including some 
fine photographs. Three papers at the 1971 National 
Heat Transfer Conference (6-8) presented velocity and 
temperature profile measurements in various configura- 
tions. One of them, Slmoneau, et al. (6), examined 
vertically Upward flow past a heated flat plate. The 
authors claimed the results Indicated the presence of 
a strong body force effect. The present paper extends 
that work to horizontal flow. 

It was the purpose of the experiment being re- 
ported herein to examine further the conclusion on the 
influence of buoyancy by comparing detailed measure- 
ments taken In both horizontal and vertical flow. 

Thus, a set of conditions from reference 6, which 
yielded a wide range of effects, was selected to be 
repeated in horizontal flow. The nominal conditions 
for this comparison were: P* » 3.74 MPa (Poo/^c “ 

1.10) s T* - 110.3 K <T./T C - 0.874); and Re x * - 
1.03xl0 6 . Data were taken at these conditions in 
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horizontal flow with the heated surface both facing 
upward and facing downward * These horizontal flow 
data are compared with the vertical flow data. 

Note should be made of a discrepancy that exists 
between this and the earlier work (6). A newer 
thermophysical property program, GASP (9), is used 
herein. It has an updated viscosity equation which 
changes the Reynolds number calculation. Reynolds 
numbers herein are about 10 to 13 percent higher. 

Any earlier data used herein has the Reynolds number 
recomputed using the newer program, GASP (9). 

EXPERIMENTAL APPARATUS 

The basic flow rig is shown in figure 1. It was 
a once-through or blow-down type system. Liqui i ni- 
trogen was transferred from a low pressure supply to 
the 0.4 m 3 high pressure dewar. The fluid tempera- 
ture was controlled by bubbling nitrogen ga> through 
the dewar. Nitrogen gas, at the desired pressure, 
was applied to the top of the fluid to force the 
fluid through the system. Syscem pressure and flow 
rate were maintained by the use of throttling valves. 
The system was heavily insulated and was precooled 
from an auxiliary supply. The system was designed so 
that the test section could be placed either in the 
horizontal or vertical position. 

The test section details are shown in figure 2. 
The test section was of rectangular cross-section of 
about 6 to 5 aspect ratio. It contained two 2.5 cm 
wide by 15 cm long Nichrome plates back to back, in- 
sulated from each other by a bakelite centerbody. 

The first centimeter was a copper electrode and was 
unheated. The primary surface was instrumented 
longitudinally with 11 Chromel-Constantan thermo- 
couples electrically isolated from the heater. Th»* 
opposite surface served as a guard heater and was 
also instrumented. Heat was generated by the dissi- 
pation of direct-current power. 

The boundary layer surveys were made with a com- 
bination pitot-static and thermocouple probe (fig. 3) 
which was traversed through the boundary layer toward 
the primary * -’ter surface. The distance from the 
primary heatek surface to the opposite wall was about 
2.4 cm, which was approximately 10 times the boundary 
layer thickness. On the back side the distance to 
the wall was 0.72 cm, about three times the boundary 
layer thickness. The temperature portion of the 
probe was an exposed ball, Chromel-Constantan thermo- 
couple encased in a sheath with an overall diameter 
of 0.25 mm. The dynamic pressure part of the probe 
was a two leg pitot-static probe of rectangular 
cross-section at the tip and had an overall tip 
height of 0.25 ran. The total pressure leg had a 
0.13 by 0.70 mm opening at the tip. The static port 
was a 0,38 mm diameter hole located six thicknesses 
from the tip to avoid separation. Measurements could 
be made to within 0.14 ram of the wall. All of the 
surveys reported herein were made 13.1 cm from the 
leading edge of the plate. 

An important feature of this design was making 
both legs of the probe identical. This was done in 
order to match the response lags in the lines, re- 
ducing the noise, and also to assure the same density 
gradient in each leg to avoid static head errors. 

The probe noise that did appear was assessed to be 
due to very small static pressure perturbations. 


Flow rates were measured with venturi flow meters. 
All system pressures were measured with strain gage 
transducers, except at the probing station where the 
system static pressure was measured with a high pre- 
cision balance type transducer accurate to ±0.05 per- 
cent. The dynamic head was measured with a variable 
reluctance type transducer. System temperatures were 
tneacured with platinum resistance thermometers. The 
heater plate thermocouples were standard Chromel- 
Constantan. All signals, Including test section volt- 
age and current, were sensed with a precision digital 
voltmeter and recorded on a digital data acquisition 
system. All measurements were checked by some kind of 
redundant Instrument. 

The system was brought up to nominal operating 
conditions and allowed to stabilize on some near- 
nominal operating condition* At this point a traverse 
of the probe from the free stream to the plate was be- 
gun. The traverse was made by moving the probe in 
about 15 discrete steps, stopping at each step to ac- 
quire data. Because of the probe noise mentioned 
earlier, at each step the signals were electronically 
Integrated for about 3 seconds. The integrated sig- 
nals and all system parameters were recorded and the 
probe was moved to the next position. Because of the 
blow down configuration with a fixed supply of liquid 
nitrogen, this procedure limited the number of points 
per traverse to about 15, which in general was suff i • 
cient to describe the profile. In addition the bulk 
fluid temperature had a tendency to drift steadily up- 
ward about 1 K during a traverse. The free stream 
conditions reported are the average over the traverse 
time. 

RESULTS 

For reference let us first examine the velocity 
distribution over the unheated plate. Figure 4 shows 
a series of velocity profiles for both gas and liquid 
flow at zero heat flux in all three geometric orienta- 
tions. In both the gas and liquid cases the profiles 
appear to be rather conventional and there appears to 
be no particular bias according to geometric orienta- 
tion. The gas has about a 1/6 power profile, whereas 
the liquid is more toward 1/5 power. In both cases, 
the edge of the boundary layer is approximately 2.0 mm 
from the heater surface. The scatter in the data from 
run to run is about 2 percent. 

Turning now to the near critics? fluid with heat 
transfer! a summary of the range of extremes that were 
obtained in these experiments is illustrated in the 
velocity profiles shown in figure 5. In the vertically 
upward flow the velocity profile exhibits a maximum 
between the wall and the free stream. On the other 
hand, in the horizontal flow with the heater facing 
downward the velocity appears to go to zero at the 
probe location nearest to the wall. The velocity 
profile with the heater facing upward 1 b more normal 
looking, however, it too has quite a steep gradient. 

With that introduction as a reference, let us 
look at a complete set of velocity and temperature 
profiles for each of the flow orientations. Figures 
6, 7, and 8 are velocity and temperature profiles for 
vertically upward flow, horizontal flow with the 
heater facing upward, and horizontal flow with the 
heater facing downward, respectively. In each case 
the nominal operating parameters are: T® * 110.3 K, 

P m - 3.744 MPa, (T c - 126.3 K, P c - 3.417 MPa), U® - 

0.660 m/sec, and Rc x ^ - 1.03xl0 6 . The surface heat 
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flu* la varied Iron 3 to IS watta/cm 2 , which In turn 
caueea the wall temperature at the probing etatlon to 
vary from ?3 K (below T c ) to over 200 K (well above 
T c ). In all eaaea the bulk of the fluid, the free 
a*-ream, la below tt. > thermodynamic critical tempera' 
rure, thua, whenever T w > T c , the critical tempera- 
ture occu.s somewhere in the boundary layer. 

In the vertically upward flow case (fig. 6), the 
effect of the boundary layer moving into the critical 
region ie most pronounced. At the lowest heat flux 
(q - 3.1) the wall temperature remains below the 
thermodynamic critical temperature and the profiles 
have a rather conventional appearance to them. When 
the heat flux is raised sufficiently to cause the 
wall temperature to exceed the critical temperature 
the velocity profiles begin to show unusual shapes, 
first gattlng quite flat and then, as heat flux is 
increased, exhibiting a maximum between the wall and 
the free stream. This effect was first observed by 
Wood (4) in a vertical tube experiment. In that case 
the presence of two walls produced two maxima and the 
profile had the shape of an M. The terminology 
Mr shaped profiles is quite descriptive and has been 
carried in the literature. The corresponding temper- 
ature profiles kept getting flatter as heat flux in- 
creased, until at the higher heat fluxes the profiles 
tended to converge on each other. These observations 
were reported in reference 6 and it was also observed 
that the significant change in behavior was accompa- 
nied by a five-fold Increase in Grashof number. 

To further examine this without the parallel 
body force the test section assembly was placed in 
the horizontal position. Looking, first, at results 
for the heater surface facing upward, the velocity 
and temperature profiles e shown in figure 7. The 
dramatic M-shaped velocity profiles of upward flow 
are not present. On the other hand, both the veloc- 
ity and temperature profiles exhibit marked inflec- 
tion points at about y » 1 mm for the higher heat 
fluxes. An examination of the temperature profiles 
reveals no special thermodynamic effect, such as 
crossing the critical temperature, occurring in this 
region. In fact, in all of the profilec, for all 
heat fluxes, the thermodynamic critical temperature 
occurs very near the wall. It is interesting to note 
that for q * 9, 12, and 15 watts /cm 2 the fluid tem- 
peratures at y - 0.14 mm are almost the same, being 
123.5, 124.0, and 125.4 K» respectively, while at the 
same time the wall temperatures are 173, 210, and 
255 K, respectively. The energy required to drive 
the boundary layer through the critical temperature 
eppeere to be quite large. This is consistent with 
the tendency of the temperature profiles to get in- 
creasingly steep near the wall and then flatten out 
away from the wall as the heat flux Increase?. Since 
buoyance is normal to the flow, the inflection of the 
velocity profiles at higher heat fluxes could be in- 
terpreted as an acceleration of the fluid near the 
wall due to thermal expansion. This question was 
raised concerning the upward flow data when It was 
first acquired. Were the so-called M-shaped profiles 
in velocity the result of buoyancy or fluid acceler- 
ation? It would appear that, while acceleration 
forces may influence the vertically upward flow, they 
are not strong enough to produce the M-shape. 

The final set of data, the velocity and temper- 
atures in horizontal flow with the plate facing 
downward (fig. 8), are alco useful in addressing the 
question of the relative roles of buoyancy and accel- 
eration. 


The most striking observation in these velocity 
profiles is that for two profiles the velocity goes to 
zero at a ^robing station a finite distance from the 
wall. Actually the transducer readings oscillated 
about zero and, since the values were well within the 
error level of the transducer, they can only be quali- 
tatively interpreted. While one could question the 
validity of saying the velocity goes to zero somewhere 
in the boundary layer, it is clear from the data 
points at slightly greater values of y that the ve- 
locity is dropping off very rapidly in the boundary 
layer at some finite distance from the wall. The tem- 
perature profiles also exhibit curious behavior. At 
q ■ 15.3 the fluid temperature sensed by the probe at 
y * 0.14 wan greater than the wall temperature. It 
appears that a highly unstable condition existed here 
and that local pockets of stagnant fluid were accumu- 
lating near the wall. This unstable condition was ob- 
served experimentally in the wall temperatures as well 
and it was very difficult to operate the facility. In 
fact, at one point the test plate over-heated to the 
point of loosening the thermocouples. Although the 
temperature profiles tend to become flatter as heat 
flux increases, the barrier near the critical tempera- 
ture is not as strong in these boundary layers. 

A couple of observations can be made by taking 
together the three sets of profiles which represent 
the three different flow geometries. First, both 
fluid acceleration and body forces play a role in 
shaping the velocity profiles, but body forces appear 
more dominant. Second, most of what is happening oc- 
curs very close to the wall, much closer than in con- 
ventional flows. One cutlous result is that at q * 

6 watts/cm 2 all three geometric orientations yield the 
same wall temperature, 145 to 147 R. Even though the 
velocity profiles are much different from case-to- 
case, the temperature distributions are quite similar. 
It is only above q ■ 6 watts /cm 2 that wide variations 
occur in the results produced by the different flow 
orientations. 

Some additional Insight can be gained by examin- 
ing the heat transfer coefficients. These are tabu- 
lated along with the system parameters for each run in 
table 1. Table I also Includes a calculated heat 
transfer coefficient. This Is based on an equation 
from Kays (10) for turbulent flow over a flat plate 
with uniform surface heat flux 

St x » 0.0307 Re x ^2 p r ;0.4 

The choice of heat transfer equation is somewhat arbi- 
trary for this discussion. What is required is an 
equation which will describe turbulent heat transfer 
in forced convection over a flat plate and the above 
is a logical choice for that. Examination of the am- 
bient gas data for horizontal flow in table I suggests 
that the Kays equation may be about 15 to 16 percent 
low in the present test application. 

The gas heat transfer coefficient for vertical 
flow listed in table I is about 32 percent above the 
average for horizontal flow. To examine whether free 
convection could account for the Increased heat trans- 
fer a standard textbook equation was selected from 
Eckert and Drake (11) 

Nu x - 0.10 (Pr.Cr Kt .) 1/3 

For the conditions of the vertical gas run this equa- 
tion yields a free convection heat transfer coeffi- 
cient, 0.00583 watts/cm 2 , which is 19 percent of the 
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forced convection value. Thus it seems possible that 
free convection augmentation could be sufficient to 
Produce to Increased heat transfer in vertical gas 


In liquid nitrogen at q - 3 watte/cm 2 , where 
the wall temperature is below critical, the results 
are qualitatively similar to the gas. In the liquid, 
however, the experimental Increase in vertical flow 
is a more modest 11 percent, whereas the computed 
free convection coefficient Is about AO percent of the 
computed forced convection value. The important point 
to note is that in all instances, regardless of geo- 
metric orientation, the heat transfer coefficients 
drop off as the heat flux Is Increased to values 
which cause the wall temperature to go above the 
critical temperature. In almost all instances these 
value b are well below those for heat transfer to sub- 
critical temperature liquid nitrogen as computed by 
standard equations such as that of Kays'; Calcula- 
tions with properties calculated on the basis of wall 
or film temperatures yield no better results. 


Another way of looking at this is to examine the 
behavior of the Grashof number. In figure 9 the 
Grashof number is plotted for the nominal operating 
conditions of the experiment as a function of wall 
temperature. The Grashof number rises steeply and in- 
creases an order of magnitude as the wall passes 
through the critical temperature region. It was 
noted in reference 6 that the unusual behavior began 
in a range 1.6 < Gr x /Re 2 < 3#0 For the present ex- 
periment this would be 1.70xlQ 12 < Gr x < 3.18xl0 12 . 
The Grashof plot shows that the body force has in- 
creased significantly in the exact range where all the 
anomalies begin. Of course, the Grashof number, as 
defined, Is only meaningful in the vertical orienta- 
tion. 


SUMMARY AND CONCLUSIONS 


An experiment has been conducted in which nitro- 
gen near its thermodynamic critical point was flowed 
past a heated flat plate. The plate was horizontal 
and data were taken for both heater facing upward and 
downward. This was compared with earlier data taken 
in vertically upward flow. The velocity and temper- 
ature profiles taken in these three orientations pro- 
vide a comprehensive and consistent set of data for 
analysts to use to develop their theories. 

In all geometric orientations the velocity pro- 
files were strongly affected when the wall tempera- 
ture exceeded the critical temperature. In upward 
flow a maximum occurred in the velocity between the 
wall and the free stream. In horizontal flow with 
the plate facing upward a distinct inflection oc- 
curred in the velocity profiles. In horizontal flow 
with the plate facing downward there were stagnant 
regions where the velocity went to zero. 

It was observed that it was very difficult to 
get the fluid in the boundary layer to cross the 
critical temperature. Thus, the temperature gradi- 
ents were very steep near the wall, much steeper than 
normal. 


Examining the profiles from the three orienta- 
tions together, it appears that both fluid accelera- 
tion due to expansion and due to body force occur. 

In the orientations where body rce can Influence 
the flow it seems to be a much more dominant force. 


Heat transfer coefficients drop off sharply as 
the wall temperature increases above the critical tem- 
perature. They are well below what one would compute 
for turbulent forced convection heat transfer to liq- 
uid nitrogen flowing past a flat plate. 

The threshold of the region of anomalous behavior 
associated with near-critical fluids begins when the 
wall temperature is slightly above the critical tem- 
pera ure. The heat flux which will produce this 
threshold temperature can be predicted by conventional 
correction. 
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TABLE I. 


Flow geometry 


HEAT TRANSFER COEFFICIENTS 


| "r Ij&c W/cm^ 


n exp» 
W/cm* K 


Horizontal - plate tip 


Gas - vertical 3.7*2 27*.*, — 

Gas - horizontal up 3.745 300.41394.: 2.75 

Gas - horizontal down 3.766 280. 5| 368.0 

a Flow appear • locally unstable. 


109.8 122.9 

110.7 143.5 

110.1 173.7 

110.6 191.3 

110.3 210.0 

0.642 

.655 

.653 

.670 

.659 

3.05 l 
6. OB l 
9. IB 1 
12.08 ; 
15.15 

L - 383 0.996 

L.395 1.032 

1 >/ 1.019 

..393 1.054 

1.392 1.033 

0.233 

.185 

.144 

.150 

.152 

0.171 

.174 

.173 

.177 

.175 

110.1 125.5 

110.4 145.2 

110.0 173.3 

110.4 210.1 

110.6 255.1 

.642 

.662 

.641 

.655 

.664 

3.07 

6.08 
9.15 

12.05 

15.22 

1.386 1.001 

1.387 1.036 
1.386 .998 
1.392 1.028 
1.395 1*046 

.199 

.175 

.145 

.121 

.105 

.171 

.175 

.171 

.174 

.176 

110.2 124.2 

110.4 147.0 

110.3 218.1 

109.7 a 200 . 1 

110.7 a 268.9 

.663 

.665 

.657 

.670 

.697 

3.07 

6.07 

9.08 
12.09 
15.29 

1.390 1.037 1 

1.389 1.041 

1.387 1.027 

1.381 1.037 

1.393 1.097 

.219 
.166 
.084 
a .133 
a . 097 

.176 

.176 

.174 

.177 

.183 

272.4 340.9 

2.61 

3.14 

.743 .928 

.0458 

.0301 

300.4 1 394.: 

2.75 

3.16 

.724 .823 

.0335 

.0291 

280.5 1 368.0 

2.74 

3.13 

.737 .931 

.0358 

.0308 
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Figure 4 - Velocity profiles for flow past an unheated flat plate. 
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Figure 9. - Variation In Grashof number over the range of the experiment for 
the nominal operating conditions. 


NASA*L«wls 



